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Abstract 

We study the connection between the appearance of a 'metastable' behavior of 
weakly chaotic orbits, characterized by a constant rate of increase of the Tsalhs 
q-entropy [J. of Stat. Phys. Vol. 52 (1988)], and the solutions of the variational 
equations of motion for the same orbits. We demonstrate that the variational equa- 
tions yield transient solutions, lasting for long time intervals, during which the 
length of deviation vectors of nearby orbits grows in time almost as a power-law. 
The associated power exponent can be simply related to the entropic exponent for 
which the q-entropy exhibits a constant rate of increase. This analysis leads to the 
definition of a new sensitive indicator distinguishing regular from weakly chaotic or- 
bits, that we call 'Average Power Law Exponent' (APLE). We compare the APLE 
with other established indicators of the literature. In particular, we give examples 
of application of the APLE in a) a thin separatrix layer of the standard map, b) 
the stickiness region around an island of stability in the same map, and c) the web 
of resonances of a 4D symplectic map. In all these cases we identify weakly chaotic 
orbits exhibiting the 'metastable' behavior associated with the Tsallis q-entropy. 
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1 Introduction 



The usefulness of the so-called 'non-extensive q-entropy' \T\ in characterizing 
the statistical mechanical properties of nonlinear dynamical systems has so 
far been demonstrated in a number of instructive examples in the literature 
(see |2] for a comprehensive review). In the present paper we focus on one 
particular property of the Tsallis q-entropy, first reported in [3|1] , and further 
explored in [5|^|[7j . These authors demonstrated that when a nonlinear dynam- 
ical system is in the regime of the so-called 'edge of chaos' the rate of increase 
of the q-entropy remains constant for a quite long time interval. Tsallis and 
the coauthors [3] argued that this behavior of the q-entropy can be connected 
to the phenomenon of a power-law rather than exponential sensitivity of the 
orbits on the initial conditions. In the case of the Feigenbaum attractor such 
a power law was observed by Grassberger and Scheunert Grassberger re- 
cently [9] questioned the meaning of the q-entropy in that particular case, but 
his arguments were convincingly rebutted in [10]. Furthermore, in the case of 
conservative systems, Baranger [6J found a constant rate of increase of the 
usual Boltzmann - Gibbs (BG) entropy in strongly chaotic systems such as a 
generalized cat map or the Chirikov standard map for high values of the non- 
linearity parameter K. However, when K is small, there is again a transient 
interval of time in which the q-entropy, rather than the Boltzmann - Gibbs 
entropy, exhibits a constant rate of increase. This phenomenon was found nu- 
merically in low-dimensional mappings when q takes a particular value (in the 
standard map q ^ 0.3 for K close to Kc = 0.97... PDIT^ . while q ~ 0.1 when 
K ~ 10~^ [13], and it was called a 'metastable state' [11]. 

The above calculations were based essentially on a 'box counting' method. 
Namely, the phase space is divided into a number of cells, and the average 
covering of these cells is found over 'many histories' [H], i.e., over large en- 
sembles of orbits, when the initial conditions are taken inside a very small 
domain (e.g. a box of size 10~^). The sensitivity to the initial conditions was 
checked by following nearby orbits with a small initial separation, (e.g. 10~^ 
or 10~^^). Finally, the value of q for which the q-entropy grows linearly in time 
was found by inspection, i.e., by trying many different values of q. 

On the other hand, the same type of 'metastable' behavior (as implied by 
Tsallis in [3]) should be recovered more simply if, instead of taking averages 
over many orbits, one calculates the time evolution of the deviation vectors 
^{t) given by solving the variational equations, together with the equations 
of motion, for the orbits inside a weakly chaotic domain. In the present pa- 
per we consider, precisely, the question of how can the 'metastable behav- 
ior'associated with a constant rate of production of the Tsallis q-entropy be 
justified theoretically by an analysis of the behavior of the variational equa- 
tions in the limit of weak chaos. This analysis is further substantiated by 
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numerical results. In fact, from the numerical point of view the method of 
variational equations is advantageous over the method of 'many histories' in 
that, when the phase space is compact, the variational equations yield the 
local rate of growth of deviations along one orbit for an arbitrarily long in- 
tegration time, while a calculation based on the integration of many nearby 
orbits reaches a saturation limit when the spreading of the orbits extends to 
the whole domain of the phase space available to them. 

The main results of our investigation are: 

a) We justify theoretically why do metastable states with a constant rate of 
increase of the q-entropy appear when a system has weakly chaotic orbits. We 
find this to be due to the time evolution of the deviation vectors, which is 
of the form ^{t) ^ at + e'^*, i.e. a combination of a linear and an exponential 
law, with a » 1 and A << 1. This law in numerical applications appears as 
producing a transient behavior for a long time interval in which the growth of 
C,{t) is almost a power law C,{t) oc t^, with p > 1. Furthermore, the exponent p 
can be associated to a q— exponent via a simple relation p = 1/(1 — g). 

b) We introduce a method by which one can calculate the q-exponent along a 
flow of chaotic nearby orbits directly from the variational equations of motion. 

c) The same method can be used as a 'chaotic indicator' distinguishing weakly 
chaotic orbits from nearby regular orbits. We call this indicator APLE (average 
power-law exponent). Its sensitivity is comparable to that of other established 
indicators in the literature such as the Stretching Numbers [H], Fast Lya- 
punov Indicator [15], Spectral Distance [16], the Mean Exponential Growth 
of Nearby Orbits [T7fT8] or the Smaller Alignment Index p3]. In section 3 we 
make a comparison of the APLE with the FLI and the MEGNO. While these 
indicators are practically equally powerful in distinguishing order from chaos, 
the APLE yields simultaneously the value of the q-exponent. 

We give examples of the numerical behavior of APLE in the domain of weak 
chaos in 2D and 4D symplectic mappings. We consider in particular: a) a 
separatrix chaotic layer in the 2D standard map, b) a stickiness domain at the 
border of an island of stability in the same map, and c) weakly chaotic orbits 
in the Arnold web of resonances in a 4D symplectic map proposed in ^U\. In 
all these cases we identify which chaotic orbits exhibit the kind of metastable 
behavior proposed by Tsallis and his associates, and which fall into the usual 
regime of the constant Kolmogorov - Sinai entropy. We also see how is this 
behavior depicted in the time evolution of the APLE, and of other indicators, 
i.e., the FLI or the MEGNO. 

The definition of APLE is given in section 2, following a theoretical analysis of 
the emergence of a transient power-law growth of deviations for weakly chaotic 
orbits in conservative systems. Section 3 presents numerical calculations in the 
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2D standard map and in the 4D map proposed in [20]. Section 4 summarizes 
the main conclusions of the present study. 



2 Theoretical considerations 

2. 1 Entropy and the growth of deviation vectors 

In order to introduce the relation of the concept of entropy to the solutions 
of the variational equations in dynamical systems, we consider the example of 
a system with a two-dimensional phase space (figured], schematic). In all the 
panels the origin 0= (0, 0) represents a hyperbolic fixed point. For simplicity, 
the eigenvectors of the monodromy matrix at O, corresponding to a mapping 
of the orbital flow after a period of time T, are assumed orthogonal, while the 
eigenvalues Ai(r), \2(T) are real and positive. The axes are parallel to the 
(orthogonal) eigenvectors. 

We consider first a non- conservative case in which both Ai(T) and A2(T) are 
greater than unity (figured^, in which we set Ai(T) = 3, A2(T) = 2). Let 
the phase space in the neighborhood of O be divided into a number of square 
cells of linear size 6. This size may represent an accuracy limit as regards our 
knowledge of the positions of the orbits due to experimental, numerical or 
other sources of uncertainty. The gray cell, of area Vq = 5^, is a set of initial 
conditions which are 'nearby' to the initial condition x(0) inside the cell. At 
a time T, the orbit and the nearby orbits are mapped to the black point and 
the six gray cells shown in figure [T]d respectively. If we ask whether the orbit 
is in the upper or lower, or, in the left or right group of cells, then after three 
questions of this type we can locate precisely in which cell the moving point is 
located, that is, we can have the maximum possible information on the position 
of the moving point in the phase space given the accuracy limit S. The number 
of questions needed to obtain this information is called Shannon's entropy 
(see [H]). Asymptotically (in the limit of a large number of cells). Shannon's 
entropy is given by the logarithm of the number of gray cells d = V{T)/V{0), 
where V{T) is the total area (or 'volume') of the gray cells at time T. This can 
be connected to the eigenvalues of the monodromy matrix at O via V{T) = 
Ai(T)A2(T)52, so that Sshan = \nd = \n{V{T)/ViO)) = InAi(T) + lnA2(T). If 
the cells are interpreted as 'micro-states'. Shannon's entropy is equivalent to 
the usual Boltzmann-Gibbs entropy Sbg = In d. The average rate of increase 
of the Boltzmann-Gibbs entropy is then given by 

Sbg Ssh an 
~Y~ ^ T ^ 
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Fig. 1. Schematic representation of the spreading of orbits in the neighborhood of 
a hyperboUc fixed point (origin) in the case of a non-conservative 2D dynamical 
system (a and b), or conservative system (c,d), see text for details. The numerical 
values in the coordinate axes are multiples of 6. 

^lnAi(r) + ^\nX,{T) = ^( ln(6(T)/6(0)) + ln(6(T)/6(0))) (1) 
where ^i(T),C^2(T) are the time T images of two initial deviation vectors 

6(0) = (5,0), 6(0) = (0,5). 

If, now, we consider the case of a conservative system, in which Ai(T) and 
A2(T) are reciprocal, the image of the set of orbits in the initial gray cell 
(figure [It, in which we set Ai(T) = 3, A2(T) = 1/3) is the dark gray area 
of the parallelogram shown in figure [T]i. However, the vertical side of this 
parallelogram is now smaller than the accuracy limit 6. Thus, if questions 
are made in order to identify in which cell the orbit is located, the search is 
restricted to a smaller number of cells d = 3. In this case we have 

5^c = lnrf = lnAi(T)^%^ = llnAi(T) = iln(6(T)/6(0)) . (2) 



5 



In the limit of large T the quantities (1/T) InAfc(T) = (1/T) ln(^fc(r)/^fc(0)), 
k = 1,2 yield the spectrum of Lyapunov characteristic exponents of the orbits 
(in the neighborhood of O). We see that in both cases, of Eq.([T]) or Eq.([2]), the 
rate of increase of the Boltzmann - Gibbs entropy Sbg/T is given by the sum 
of positive Lyapunov characteristic exponents of the orbits, which according 
to Pesin's theorem |22], is equal to the Kolmogorov-Sinai entropy of the orbital 
flow in the neighborhood of O. This means that the Kolmogorov- Sinai entropy 
is, in fact, a measure of the rate of change of the entropy rather than a measure 
of the entropy itself. 

The following is a more precise treatment of the previous schematic picture. 
Consider a partitioning of the n-dimensional phase space 7W of a coservative 
system into a large number of volume elements of size 5" for some small S > 0. 
Let x(0) be the initial condition of an orbit located in a particular volume 
element Vq = ^oi^o2---^on, where ^ofc, k = 1, ...,n are the linear dimensions of Vq 
in a locally orthogonal set of coordinates in the neighborhood of x{0). Without 
loss of generality, we set all the initial values ^ofc equal, i.e., S = Vq^^ = ^ok, 
\/k. All the orbits with initial conditions within Vq are called 'nearby' to the 
orbit x{t) with initial condition x(0). Because of the volume preservation, the 
orbital flow defines a mapping Vq V{t) of the volume Vq to an equal volume 
V{t) at the time t. We want to find an estimate of the covering of the cells 
of A4 by the volume V{t) in terms of the variational equations of motion. To 
this end, let 

m = Dtio (3) 



be the solution of the variational equations for an initial deviation vector 
acted upon by a linear evolution operator Dt determined solely by the orbit 
x{t). Let Ck{t) be the images of ^ofc, k = l,2,...,n under the action of the 
operator Dt- The vectors {^^(t)} form a complete basis of the tangent space 
to M at the point x(t) iff {^ofc} form a complete basis of the tangent space to 
Ai at the point x{0) and rank{Dt) = n. It is possible to obtain an orthogonal 
basis {C,k{t)} starting from {^^(t)} via the Gramm-Shmidt procedure [23j. The 
new basis is obtained by the recursive relation: 



m=m-j:i^kit)-ut))j%l, . (4) 

The volume V{t) is given by V{t) = ^i(t)C^2(t)---(,nit). We reorder this basis 
by decreasing length of the vectors C,k(t). Let Vc{t) be a coarse-grained volume 
equal to the total volume of all the cells visited by the orbits in V{t). Vc{t) is 
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determined by only those vectors with lengths greater or equal to S, that is 



Ve{t)=^,it)Ut)-Ut)S''-"' (5) 

where m is defined by the condition m = sup{m' : ^k{t) > S for all k = 
1, . . . , m'}. It follows that m = n when t = and 14(0) = 1^(0) = 5". 

The Boltzmann - Gibbs entropy is defined by 

SBG{t)^\nW{t) = \n^^ (6) 



where 



yields the number of cells (or 'micro-states') occupied by Vc{t). 

The rate of growth of the entropy ([6]) for a set of nearby orbits is related to 
the spectrum of the Lyapunov characteristic exponents of the reference orbit 
x(t). The Lyapunov characteristic exponents are given by 

A,^limiln^=limlln^ (8) 



ioT k = 1, . . . ,n, provided that the limits exist. Let Ai, A2, ...Am, < m < n/2 
be the set of positive exponents of the spectrum ([H]). According to Pesin's 
(1978) theorem, their sum is equal to the Kolmogorov- Sinai entropy Sks of 
the fiow of orbits 'nearby' to x{t) [2H[25] . i.e. 

5^5 = Ai + A2 + ... + Am . (9) 



The average rate of increase of the Boltzmann-Gibbs entropy up to the time 
t is given by SBG(t)/t- In view of Eq.Q, the limit of this rate, for t — > 00 is 

o 1 m 

lim ^ = lim - In = ^ Afc = Sks (10) 



Thus, in the limit t ^ 00 the Kolomogorov-Sinai entropy is equal to the 
asymptotic value of the mean growth rate of the Boltzmann-Gibbs entropy. 
Precise numerical examples of this relation were given in the case of low 
dimensional mappings pBfT] . However, the orbits may exhibit transient or 
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'metastable' states for long time intervals before reaching the limit (fTOl) . Such 
states are characterized by a constant growth rate of the Tsallis q-entropy, to 
which we now turn our attention. 



2.2 Tsallis entropy and the Average Power Law Exponent 



The time evolution of the q-entropy [T] for an ensemble of orbits with initial 
conditions within the volume Vq is given by 

s,i,^'mg-i. (11) 



In this equation g is a constant parameter, known as the 'entropic index' q. 
Dividing Sq by t/ti, where ti is a transient initial time of evolution of the 
orbits, and substituting W from ([7]), yields the mean rate of evolution of Sg 

% = 7T7r^. • (12) 



For every C,k, ^ = 1, 2, ...m, we define an Average Power Law Exponent (APLE) 
Pk according to 

Ut)=Ck{ti)i^r', k = l,2,...m . (13) 
tl 



All the Pk are, in general, functions of the time t, and the value of Pk{t) 
yields the average logarithmic slope (or power-law exponent) of the evolution 
of ^ (t) in the time interval from the time ti up to the time t. Furthermore, in 
conservative systems we have pi + P2 + ■■■ + Pm > 0, since, by the preservation 
of volumes, the components ^k{t) cannot be all decreasing functions of the 
time. 

In view of the definition of the APLEs 0131) . equation (I12p ) takes the form 

^1 ^ ((^ \(Pl+P2 + -Pm)(l-g) _ 1\ (1A\ 

t/t, {t/t,){l-q)^H,^ ' ■ ^ ^ 



In the limit t — > oo the quantity ^ tends to a non-zero finite value only if a) 
the PiS take constant limiting values, and b) the entropic index q satisfies the 
relation (pi +p2 + ...pm)(l ~ ?) = 1- In all other cases, ^ tends either to zero 
or to infinity. If the deviations ^fc(^) grow asymptotically as a power law, then 
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condition (a) is satisfied and the mean rate of increase of the Tsalhs entropy 

S 

jjj^ tends to the sum of the positive APLEs 
S 

hm —J- =p^+p2 + ...pm (15) 

t^oo t/ti 



for the value of q given by 

q = l ^ . (16) 

P1+P2 + ■■■Pm 



In that case, if pi is by definition the maximum of all the APLEs, this exponent 

s 

can be used lower bound of the limit of 777-, i.e. 

t/ti ' 

lim 4^ > pi . (17) 



In practice we can use Eqs. ffTSl) . f|T6|) . or f|T7j) for a long but finite time t in order 
to estimate the average value of the q-exponent in the interval from ti and t, 
provided that this value is almost constant in this interval. Furthermore, the 
ratio of the length of the deviation vector C,^(t) = J2T=i ^k(J^) ^ with respect 
to the length of this vector at ti can be evaluated from the equation 



e{t) ^ ait) ^ ,2, t 



\2p 



k 



where equation f|T3l) has been used and 



,2 



with 



k=l 



Equation ( fT8|) can also be written as 

^ = ilf'-i^i + tfkilr'^'^'^'] (21) 
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Since pi — pk is positive for 2 < k < m, the sum inside the square brackets 
in the last expression tends asymptotically to zero for t » ti. Thus we can 
write 



APLE = p 




Pi + 



21nf 



(22) 



21nf 



i.e in the limit of t oo APLE tends to pi. 

In two dimensional maps or in the Poincare surface of section of 2D Hamilto- 
nian systems we have m = 1, i.e. there is only one positive exponent pi which 
is the limit of APLE p. In these cases the APLE p is also the limit of the mean 
rate of growth of Tsallis entropy 5*^ with q = 1 — 1/p according to (|T6l) . 

2.3 The time evolution of APLE for weakly chaotic orbits and the appearance 
of 'metastable ' states 

In the sequel we are interested in the behavior of the deviations ^(t) for or- 
bits in the border of a single resonance domain of a nonlinear Hamiltonian 
system of n degrees of freedom. In the integrable approximation this border is 
separatrix-like. The Hamiltonian in resonant normal form (see e.g. |27j) reads: 



in action -angle variables ( J^, J2, . . . , Jn) and (V', 4>2, ■ ■ ■ 4>n)- The resonant vari- 
ables (J^, '?/'), E 71, ip & obey a pendulum dynamics. The point 
{J^,ip) = 71 defines a foliation of simply hyperbolic n — 1 dimensional invari- 
ant tori of fl23l) labelled by the constant values of Jj, i = 2, . . . n. In particular, 
when n = 1 there is only a 0-dimensional torus, i.e. an unstable equilibrium 
point, while if n = 2 the tori are 1-dimensional, i.e., a family of unstable peri- 
odic orbits. The remaining phase space is foliated by n-dimensional tori. The 
deviations ^(t) on these tori grow in general linearly ^{t) ^ ^0 + A t, where 
A is a measure of the frequency differences between orbits on nearby tori. We 
may assume that the derivatives \duJi/dJj\ = IO'^Hq/ dJidJj\ are bounded from 
above for all i,j = 2, . . . ,n. This, however cannot be true for the frequency 
associated with either librations or rotations in the plane (J^, ■?/'), since the 
derivative of either the libration or rotation frequency with respect to the res- 
onant action tends to infinity when the action tends to its limiting value on 
the separatrix, i.e.: 




(23) 




(24) 
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where, e.g., for librations 

l/'maa;(-E) 



ME) = 2 J J4E,ij)di: (25) 



is the hbration action of an orbit labeled by the energy E = — Uq cosip 
and iprniniE) , ipmax{E) are the limiting values of ip in the domain in which the 
latter equation has solutions. Furthermore Ur{Jr) = dE/dJr, and Jr,separatrix = 
Jr{E = uJq). a similar expression is found in the case of rotations, but with 
different limits of the integral (1251) . Near the separatrix, the value of A in the 
linear term of the growth of deviations is determined essentially by the value of 
the derivative a = \du!r{Jr)/dJr\, which is finite, but large. Thus the growth is 
essentially determined by the growth of the projection of the deviation vector 
on the plane {Jr,4>r), where 0^ is the angle conjugate to Jr according to the 
previous definitions. The variational equations for ^ = (AJ^, A0r) read 

AJr = 0, A(f)r = ^^Jr (26) 



yielding the solution AJ^(t) = C = const., |A0r(t)| = |A0^(O) + C a t|. If 
C = the deviations remain constant ^(t) = |A0r(O)|, while if C 7^ they 
grow linearly in time. In numerical applications we usually select a random 
orientation k. = A0r(O)/AJr(O) = A0r(O)/C for .^o, so that, in general C = 
Co{l + ^2)^1/2 is of order C = 0(^o) 7^ and ^(t) ^ ^o(l + at). The time 
behavior of the APLE for a deviation growing as ^(t) ~ ^o{l + at) is shown 
in figure [2^, when a < 1 (curve (1)), or a > 1 (curve (2)). In fact, when t is 
large and C = 0(^o) we have ^{t) ^ ^0 a t, thus 

p^m!m = ^^l^a + lnt) . (27) 
int mt 



If a < 1 we have In a < and p tends to the value p = 1 from below. On the 
other hand, if a > 1 we have In a > and p tends to the value p = 1 from above. 
Thus, far from the separatrix the time evolution of p is like in curve (1) of 
figure [2^, while close to the separatrix it is like in curve (2) of the same figure. 
Numerical examples of this behavior are given in section 3 below. Finally, 
independently of the distance from the separatrix, if the initial vector ^0 is 
almost tangent to an invariant curve, we have C << ^q. In that case we have 
^ (t) ^ ^0 S'lid p ^ for t < 1/C, while for t > l/C p increases approaching 
asymptotically the value 1 from below, independently of the value of a (figure 
12^, curve 3). Numerically, we find that this behavior can only happen when 
the angle between ^0 and the tangent to the invariant curve is small (below 
two degrees, see section 3). Near this value there is a continuous transition in 
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Fig. 2. The time evolution of p vs. log^Q t for (a) regular orbits, and (b) weakly 
chaotic orbits, according to the simplified model of Eq. (|3ip . In (a) we have (1) 
i = v^A J2 + A,^2 and = 1, a = 0.5, e ^ 0, (2) = 1, a = 2, e ^ 0, and (3) 
^ = 1 + Y^A J2 + A(^2 ^ith Co = 10"^, a = 2, e ^ 0. In (b) we have (1) ^ as in 
(a) with Co = li = 0.5, e = 10~^, (the Lyapunov exponent and upper crossover 
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a narrow interval of values of cf) from the curve (3) to the curve (2). Practically, 
when the initial orientation of Co is selected randomly, in the great majority 
of cases we encounter for regular orbits only the cases (1) and (2) of figure [2^. 

We now examine the time evolution of the APLE in the case in which a 
Hamiltonian perturbation is introduced, namely 



H- 



■ Hr{Jr) + Hq{J2, . . . , Jn 



ei^i(J,0;e) . 



(2J 



In Eg. 0281) we assume that an optimal resonant Birkhoff normal form has al- 
ready been constructed (e.g. This means that the action - angle variables 
in fl28|) are obtained through a near-identity transformation from the original 
action - angle variables of the unperturbed Hamiltonian. Furthermore, in the 
Nekhoroshev regime the size of the perturbation e in Eq. fl28l) is exponentially 
small in the quantity Ignoring the small components of the deviation 

vector normal to the resonant plane ( J^, (pr), the new variational equations of 
motion read: 



d{AJr 
dT 



-e 



12 



^ - e(^A6 + ^ A J ) + ^AJ (29) 



While a detailed exploration of the solutions of Eqs. fl29l) can only be made 
after Hq and Hi are known, we can explore the basic behavior of the solutions 
close to the separatrix limit by the following heuristic analysis. Close to the 
separatrix we have \d^Hr/dJ^\ = \du!r/dJr\ » 1, while dur/dJr < 0. As- 
suming that all the partial derivatives of Hi in (!29l) have 0(1) average values 
over the basic periods of motion, we introduce the following simplified model 
yielding essentially the behavior of the resonant components of the deviation 
vector: 



— — — = -e{AJr + Acpr 

d{A(j)r 



dt 



e 



{AJr + A(pr) - aAJr (30) 



for a > 1 and e << 1. For e = 0, Eqs.( l30i) take the form of the equations 
(l26i) of the integrable case. If we choose initial conditions perpendicular to the 
invariant curves of Hr, i.e., AJr(O) = C = A0r(O) = 0, the solution of (!30!) 
reads: 



AJr{t)=U cos\i{{ae)'/h) - (-)'/' sinh((ae)^/2^) 



a 

AMt)=^o[ilf'-i-f')smhiiaey/H) . (31) 

The asymptotic analysis of (13T|) yields now the reason why we observe a 
'metastable' behavior in the time evolution of deviations. If e << 1, then 
for times t < (ae)~^/^ we can consider both e and (ae)^/^t as small quanti- 
ties. Then, the first of equations ([21]) yields an almost constant term AJr ~ 
^o + 0(et), while the second equation yields a linear behavior | A0r| ~ ^o\^ — ci\t- 
This is similar to the case of regular orbits. We thus have a behavior simi- 
lar to the curve (2) of figure [2^. However, when t > (ae)^^/^ the exponential 
behavior becomes dominant AJ^ ~ A0 ~ exp((ae)^/^t), with Lyapunov expo- 
nent A = (ae)^/^. This behavior is exemplified in figure [2)d, in which we plot 
^{t) = (Aj2(t)+A02(t))^/^ with AJr{t), A(j)r{t) given by Eq.dSiD, for different 
values of a and e. We see that the combination of the linear and exponential 
laws creates a 'plateau' of nearly constant value of p between an initial time 
tic ~ 10 and a second time which is essentialy given by t2c ~ A^^. These 
times are called 'crossover times'. In the interval tic < t < ^2c the deviation 
vector ^(t) grows almost as a power law ^(t) oc t^ for a nearly constant value 
of p > 1. We stress that the real time evolution ^(t) oc at -f exp(At) is mathe- 
matically different from a power law, and only a numerical resemblance to a 
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q-exponential is systematically obtained for specific time intervals. As shown 
in figure [2]d, the duration of the 'metastable' behavior decreases when A in- 
creases, while the value at which p is stabilized in the interval tic < t < ^2c 
increases as a increases. An analytical estimate of the plateau is obtained by 
noticing that the leading terms of Eq. fl3T]) (for small and large times) yield 
the time evolution of ^{t) more precisely as ^(t) ~ 1 + at/2 + (a/2A) (^e^^ — 1^ . 
Taking p ~ In ^/ In t we readily find the profile of the function p{t) in the 
neighborhood of a characteristic time to = given by the following bounds 
for the first and second derivatives of p (for A -C 1): 

dp ^ ^/2 + 1/e ^ ln(A + ae/2) 



dt - \ |lnA| (InA)^ 



fp^ 2 ( 1-2/e-l/e' 1 + ln(a/2) - 2/e 21n(ae/2) \ 

dt^ ~ \ |lnA| (lnA)2 ^ |lnA|3 J " 

The second of the above equations implies that the function p{t) is convex at 
the time to = 1/X, so that the variations of p over intervals At around to are 
bounded by the 0(A/| lnA|) estimate for the first derivative, namely: 



For example, if A = 10~^, a 1% variation of the value of the APLE Ap = 0.01 
can only occur in an interval At ~ 0.01| ln/A|/A, or At ~ 10^. This value 
marks the extend of the plateau, which is a considerable fraction of the time 
to ~ 1/A. The value of p in this plateau is estimated as: 

ln(ae/2) 

Pplateau ^ 1 + (33) 



The estimate fl33p is found to be in good agreement with the numerical values 
of p given, e.g., in the examples below. Recalling that a is essentially a measure 
of the derivative of the frequency \dur/dJr\, we expect that, in a chaotic layer 
resembling a separatrix, but with some thickness, p increases as we approach 
to the center of the layer, since a increases abruptly, while p is smaller near the 
edges of the layer. This is precisely what is found in numerical experiments, 
as analyzed in the following section. In particular, if one defines an average 
value < p > over the whole chaotic layer, this yields an average value of the 
entropic index q = 1 — l/< p > corresponding to the metastable behavior of 
the orbits in this chaotic layer. 
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3 Numerical applications of the APLE 

3.1 Separatrix layer in the 2D Standard Map 

In the sequel we consider examples of numerical calculations of the APLE = p 
in discrete conservative systems. The time t in Eq. fl22p obtains discrete values 
t = 1, 2, .... In order to avoid a singular value of p when t = 1, we set ti = 2 
and we calculate p for t > ti, with p = when t = ti. 

We are particularly interested in the time evolution of p for orbits near or 
within a domain of weak chaos. A basic example is provided by the 2D stan- 
dard map \28\ : 



Un+l = yn + TT sin i^TTXn), 

Xn+l=Xn + yn+l, (34) 

where x„, ?/„ are given modulo(l) in the intervals [0,1) and [-0.5,0.5) respec- 
tively. Figures [3^,d show the phase portrait of the map fl34p for K = 10~^ 
(figure [3^) and K = 10"^ (figure [311). The thin solid lines in these figures are 
the unstable asymptotic manifolds emanating from the unstable periodic orbit 
Pu = (0, 0). The manifolds are calculated by taking many initial conditions in 
a small segment of length 10"*^ along the unstable eigendirection given by the 
monodromy matrix at Pu. Both phase portraits show the basic structure of 
the standard map for small K, i.e., a thin separatrix layer of weak chaos that 
separates a librational from a rotational domain. When K is small the latter 
domains are filled almost entirely by invariant curves. 

Figures [3]D,e show the typical time evolution of APLE for regular or chaotic 
orbits. The dotted curves give the time evolution of the APLE for a regular 
orbit inside the domain of librations (initial conditions: (xo,?/o) = (0.0,0.2), 
close to the stable periodic orbit at (0,0.5)), for K = 10^^ and K = 10^^. The 
initial deviation vector ^(0) in both figures is chosen to be nearly perpendicu- 
lar to the invariant curve passing through {xo,yo). The temporal behavior of 
APLE in the dotted curves of figuresj3)D,e is typical of regular orbits, i.e., the 
APLE grows slowly tending asymptotically to p = 1 from below (as in curve 
(1) of figure [2^). The oscillations of p around its local mean value are due 
to oscillatory variations of the component of the deviation vector ^(t) locally 
tangent to the invariant curve. In fact, if we approximate the invariant curve 
by an ellipse, it can be shown that the amplitude of the oscillations of p is 
proportional to the axial ratio of the ellipse. Furthermore, a plot of the finite 
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Fig. 3. (a) phase portrait of the standard map for K = 10^^ (b) The time evolution 
of APLE for a regular orbit in the libration region of (a), and for a chaotic orbit 
in the thin separatrix layer, with initial conditions along the unstable manifold of 
the periodic orbit (0, 0). (c) the time evolution of the finite time Lyapunov number 
\ogiQXif) foi^ the same orbits. (d,e,f) same as in (a,b,c), but for K = 10^^. 

time Lyapunov number 

X(t) = -ln|^| (35) 



for the same orbits (figures [3b, f, dotted lines) shows also the behavior expected 
for regular orbits, i.e., x(^) falls asymptotically as for large t. 

Now, the thin solid curves in figure [3]D,e show the behavior of APLE for chaotic 
orbits inside the separatrix layers of figures [3K,d. In this case we take the 
initial conditions on the unstable manifolds of P{/, a fact ensuring that all 
the consequents of the chaotic orbits are on the same manifolds. The initial 
deviation vector ^ (0) is chosen perpendicular to the unstable manifold. We can 
immediately notice the difference in the time behavior of APLE for these two 
orbits. In the case of the solid curve of figure [3^ {K = 10~^), the APLE grows 
initially crossing the value p = 1 at a short crossover time tic = 25. However, 
after this crossing the APLE describes a number of oscillations around a mean 
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value that remains systematically above unity, up to a second crossover time 
= 38.000. We find < p >= 1.9 in the time interval tic < t < t2c- As shown 
in figure [3]F, the value at which the finite time Lyapunov number X{t) stabilizes 
is LCN = A = 10~^'^. We thus see that the crossover time t2c is essentially 
given by t2c ~ A^^. On the contrary, in the case of the solid curve of figure [3)d 
{K = 10"^), the APLE grows from the start indefinitely, as expected for an 
exponential growth of deviation vectors, and there is no visible 'metastable' 
behavior in the time evolution of p. In that case the Lyapunov number (figure 
[3t) is rather large A = 10~^, and the corresponding crossover time t2c ~ 10^ 
is of the same order as tic, i-e., extremely short to produce any visible effect. 

As in the case of regular orbits, the oscillations of APLE around a local mean 
value in figures [3)D,e are due to the oscillatory behavior of the component of 
the deviation vector ^{t) which is tangent to a theoretical separatrix passing 
through the center of the separatrix chaotic layer. In this case we find that 
the first maximum value of p in figures [3]D,e occurs when the orbits pass close 
to the first homoclinic point of the unstable manifold emanating from Pu and 
the stable manifold emanating from the image of Pu modulo 1. 

In order to check the dependence of the time evolution of the APLE on the 
initial orientation of the deviation vector ^(0), corresponding to the value of 
the constant C in the theoretical analysis of subsection (2.3), the following nu- 
merical test is performed: Starting from an arbitrary initial orientation of the 
deviation vector .^(0), the orbit and the variational equations are integrated 
for a long time t = r. It is well known (see, for example [29]) that the deviation 
vector ^(t) evolves so that it tends to become tangent to the invariant curve, 
in the case of a regular orbit, or parallel to the direction of a nearby unsta- 
ble asymptotic curve (invariant manifold), in the case of a chaotic orbit. Let 
{xr, Ut) be the position of the orbit at t = r and {dxr, dyr) be the components 
of the deviation vector at this time. Taking r long enough so as to ensure that 
the limit of tangency was reached down to the numerical precision level, we 
use [xr, Ur) and {dxr, dyr) as the initial conditions of a second orbit and of its 
variational equations. We then compare the time evolution of the APLE for 
this deviation vector and for a second deviation vector associated to the same 
orbit, but with orientation forming an angle with {dx^, dyr)- 

Figure m shows examples of the dependence of APLE on the initial orientation 
using three different orientations, namely ~ 0, = 2° and (j) = 90°, in the 
case of a regular orbit with A' = 0.1 (first row in figure H]), and in the case of 
a chaotic orbit with K = 0.1 (second row in figure H]). Clearly, the evolution 
of APLE is sensitive on (p only for small values of (p. For example, in the case 
of the regular orbit, when ~ (figure H^) we have A J = C ~ and the 
corresponding term in the solution of Eq. (!26|) is suppressed. Thus ^ (t) simply 
makes oscillations around the mean value < ^(t) >= ^(0), and the APLE 
remains close to a zero value even after t = t = 10^ iterations. On the other 
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Fig. 4. The time evolution oip along a regular orbit with initial conditions (xq, yo) = 
(0.5000046134885427,0.08097827536667734), for if = 0.1 and three different devi- 
ation vectors forming an initial angle (a) 0*^, (b) 2^^, (c) 90^ with the tangent to 
the invariant curve. The corresponding evolution for a chaotic orbit (initial condi- 
tion (0.5029208124108719,0.09994913824229859) is shown in (d-f), for K = 0.1. The 
bold dashed line in (d) is the theoretical solution corresponding to an exponential 
law C{t)/Co = e^* with A = 10-^•^ while in (e),(f) we superimpose to this a linear 
term at with a = 2. 



hand, when is equal to only = 2°, the time evolution of APLE becomes 
already very similar to its typical behavior, as concluded by a comparison with 
the case (p = 90°, corresponding to an initial deviation vector perpendicular 
to the invariant curve. The same phenomena apply to the case of weakly 
chaotic orbits as in the second row of figure HI In that case, the asymptotic 
limit for all three values of the initial angle is an exponential growth of the 
corresponding deviation vectors, leading to a constant limit of the Lyapunov 
number A ~ 10~^'^. When = the 'metastable' behavior does not show up 
in the time evolution of the APLE. However, this behavior is clearly seen when 
> 2°, and it lasts up to the crossover time t2c ~ 10'' > 10^'^ = 1/A. In fact, 
plotting the theoretical solution (13T1) corresponding to a choice = 90° in this 
case shows a good agreement with the numerical results for both = 2° or 
= 90°. 
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3.2 The APLE as a chaotic indicator measuring the entropic q-index. Com- 
parison with FLI and MEGNO 



In theory, the use of APLE as a 'chaotic indicator' distinguishing regular from 
chaotic orbits is straightforward. In the case of regular orbits the value of 
APLE tends to p = 1 as t — oo, while in the case of chaotic orbits we have 
p — *■ cxo as t — ^ oo. In practice, however, one can only evaluate p over a finite 
integration time T and a numerical indicator of chaos is considered as efficient 
if this time is small. An example of efficient chaotic indicator that is widely 
used in the literature is the Fast Lyapunov Indicator (FLI) [13] in its revised 
form [5U]. Let ^(t) be the length of the deviation vector of an orbit. The revised 
form of FLI reads: 

FLI = sup{ln {m/mi < t < T} . (36) 



In the case of regular orbits the deviations grow linearly, so that In (^ (t) /^(O)) ~ 
Int. We can thus set a threshold value, say FLIq{T) = In 10+lnT, correspond- 
ing to a deviation vector larger by a factor ten from the one corresponding to 
the linear growth of deviations up to the time t = T. Then, if FLI > FLIq the 
orbit is called chaotic, otherwise it is called regular. In fact, if we are close to 
the border of a separatrix chaotic layer, we have ^{t) = ^o(l + o, t) for regular 
orbits, with a » 1 (subsection 2.3). Thus In (^(t)/^(0)) ~ Ina + Int. If we 
choose the threshold of FLIq as above, then, if a > 10, a regular orbit can 
be erroneously characterized as chaotic. Thus, an improved formula for the 
threshold value is 

FLIo{T) = In 10 + In a + In T 

where the value of a can be estimated by the value of ^(tic) for a time 
1 << tic « T, since, if a » 1, we have ^(tic) — ^o^tic- Then, the or- 
bit is considered as chaotic if FLI > FLIq at the time T. 

Suppose now that the orbit is weakly chaotic, so that a 'metastable' behavior 
exists for which ^{t) ~ ^(tic)(V^ic)^, with p > 1 (figure [5^, in which a plateau 
is formed at about p = 1.4, from tic = 10^ to t2c = 10^). As shown in figure [5)d, 
the value of In (^(t)/^ (tic)) crosses the line of FL/o(t) (solid straight line) at 
about the same time (T = 10^) when the slope of the quantity In (^(t)/^(0)) 
vs. Int crosses the value p = 1.4 upwards, i.e. t2c ~ T. This means that the 
APLE yields the characterization of the orbit as chaotic at the time t2c which 
is of the same order as the minimum time T needed by the FLI. In fact, after 
t > T the exponential growth of deviations becomes dominant and the slope 
of In (C(t)/^(0)) vs. Int tends very quickly to infinity. 
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Fig. 5. (a) A theoretical example of evolution of the APLE when an oscillatory 
term is added to the solution ([3T]) . namely ^ = -\/AJ^ + A(^(l + cosO.lt) with 
A J,., A(/),. given by Eq.jH]) and = 1, « = 30, e = 10^^°. (b) The evolution 
loSio(|(^) l°§io^ same example. The dashed line corresponds to a linear 

law In^(t) = Ina + Int + In 10, while the continuous line indicates a power-law 
^(t) oc with (p = 1.4, i.e., as given by the approximate plateau of (a). 




Fig. 6. (a) The maximum value Pmax in a time interval 2 < t < 10^ as a function 
of the initial condition yo of the orbits in a segment along the line xq = 0.5 crossing 
the separatrix domain of figure [3ti. The peak value is for the orbit with initial 
conditions (xo,yo) = (0.5,0.0998071243), and initial deviation vector = {-^, :^)- 
The evolution of the APLE and of log^^g x{t) foi" the same orbit are shown as solid 
black curves in (b) and (c) respectively. The dotted black and the gray curves in 
the same plots correspond to orbits with initial conditions (0.5, 0.099807124) and 
(0.5, 0.09980712) respectively 



In practice, we found that the location of thin chaotic layers in resonances can 
be determined using APLE in a way analogous to the Eq. ( l36l) . namely: 

Pmax = sup [p = l^^jTJ^^ ' tl<t<T} . (37) 
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Fig. 7. (a) The time evolution of logio(|^), ^^'^ i^) time evolution of the 
MEGNO (y) for the same orbits (gray and black dotted) as in figure [6)3,0. 

In order to probe numerically the sensitivity of Pmax to thin chaotic layers, 
figure shows the variation of the value of Pmax, for 2 < t < T = 10^, along a 
segment of the line of initial conditions Xq = 0.5 passing through the separatrix 
chaotic layer of figure [3^. As the center of the resonance is approached, the 
value of Pmax increases abruptly, the peak value marking clearly the center 
of the chaotic layer. The peak value of p corresponds to an orbit with initial 
conditions {xq, Uq) = (0.5, 0.0998071243) (the initial deviation vector is taken 
as ^0 = {^-i Tf))" Si^^^ chaotic layer is very thin, the behavior of APLE 
in a thin domain including this orbit is very sensitive on the choice of initial 
conditions. Thus, the central orbit has Lyapunov number LCN ~ 10~^ (figure 
int), and it shows only a small plateau in the values of p at p ~ 3.5, for about 
900 periods, i.e., from t = 10^ to t = 10^ (figurelHb)- However, if we only change 
the initial value z/q by cutting the last digit, {xo,yo) = (0.5,0.099807124), the 
new orbit has a much smaller Lyapunov number {LCN ~ 10~^, and the APLE 
forms a conspicuous plateau at p ~ 1.5 which lasts for 10^ periods. Finally, 
if we cut one more digit in i/q, (xo,?/o) = (0.5,0.09980712), the orbit shows 
no sign of chaos up to t = 10^, although the plateau formed in the values of 
APLE (gray curve in figure Eb) is definitely at values p > 1, an indication that 
the orbit may finally be chaotic. We notice that these results depend also on 
the machine precision, and different runs with different machines or precision 
levels yield qualitatively the same picture as in figure O but for somewhat 
different choice of initial conditions. 

Figure [7^ shows, for comparison, the identification of the same orbits as in 
figure [6] by the FLL Clearly, the two methods have a similar sensitivity, i.e., 
T ~ tsc = 10^ for the dotted black curve of figure [7^. Figure [Tb shows the 
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behavior, for the same orbits, of yet a different chaotic indicator: the MEGNO 
[T7] . The definition of MEGNO for continuous systems is: 



whereas, in the case of mappings, ^(t) is replaced by the finite difference of 
the deviations ^ at successive time steps. If the deviations have an average 
power law behavior ^(t) ~ V in the interval < t < T, Eq. fl38l) yields the 
value of MEGNO Y = 2p. However, this is not so when the power-law is 
transient and approximate, i.e., it comes from a combination ^{t) ^ a t + e^*. 
In fact, for regular orbits we readily find Y{T) = 2(1 — ln(l + aT)/{aT), thus 
F(T) < 2 for T > 1/a even if a >> 1, that is the MEGNO cannot yield a 
power-law exponent p > 1. Furthermore, we find that the numerical behavior 
of the MEGNO shows variations which are about one order of magnitude 
larger than those of APLE in the whole interval of 'metastable' behavior of 
the orbits. This is exemplified in figure Uh, referring to the same orbits as in 
figure [6] or [7^. The fast oscillatory variation of the MEGNO, of amplitude 
Ay ^ 10, shown in figure [Tjo results in that, while the maximum value of Y 
is positive and above the threshold Y = 2, the minimum value is negative. 
The mean value < y > in the transient interval of time is below 2. On the 
other hand, not only the corresponding values of the APLE (figure |6]d) are 
always positive, but they are clearly above the threshold p = 1 during the 
whole interval < t < T. Thus, the identification of the metastable behavior 
is much more clear by the APLE than by the MEGNO, and the APLE can 
be used to obtain a useful quantitative estimate of < p >, and thus of the 
entropic index q = 1 — 1/ < p >. This is shown in figure ([8]) , which yields the 
average value of < p > along an orbit in the time interval < t < T = 10^, 
as a function of the orbit's initial condition, for all the orbits in the same 
scanning of a thin chaotic layer as in figure [61 The general structure of this 
diagram is as in figure [H^. However, only a subset of these orbits have < p > 
larger than unity. These orbits yield an average value of the values of < p > 
equal to < p > = 1.07, corresponding to an entropic index q = 0.065. 

3.3 Stickiness region 

A case of particular interest, in which weak chaos emerges, is in the stickiness 
region at the border of an island of stability. The stickiness of the orbits is 
due to the existence of one or more cantori (see [31] for a review). The cantori 
arise from the destruction of KAM tori, which, according to Greene's criterion 
[32] happens when the stable periodic orbits with rotation numbers forming 
sequences with limit equal to the rotation number of a torus become unstable 
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Fig. 9. The time evolution of the APLE for some orbits belonging to the sticky 
region of an island of stability in the standard map for K = 5 (see [37j). The 
initial conditions of the orbits are (a) (xo,yo) = (0.64337,0.36), (b) (xo,yo) = 
(0.6433684,0.36) and (c) (xo,yo) = (0.643363,0.36). 

j33|l34|l35|36] . Since the dynamics at this hmit is very close to hyperbohc, we 
expect that weakly chaotic orbits in the stickiness domain exhibit a transient 
metastable behavior for times of the order of the stickiness time. 

In order to study this phenomenon, we take initial conditions at the border of 
an island of stability as in [37J, namely we consider initial conditions in a line 
segment with i/q = 0.36 and 0.64324 < Xq < 0.6434 (as in figure 7 of [37]). The 
stickiness time is particularly high (it can be larger than 10^ periods) when 
Xq > 0.64336. By integrating many chaotic orbits in this domain, we found 
three kinds of different behavior in the time evolution of the APLE, shown in 
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figures ([9])a,b,c respectively: (a) An orbit may form no 'plateau' beyond p = 1 
until the orbit escapes in the large chaotic sea outside the stickiness zone, (b) 
An orbit may form one plateau lasting for times of the order of its stickiness 
time, e.g. the orbit of figure [Hb has a clear metastable behavior, with p ~ 1.3 
for 7000 < t < 25000, while the orbit escapes after about 10^ periods, (c) An 
orbit forms more than one 'plateaus' before escaping (the orbit of figure [Qt 
forms two large plateaus at p ~ 1.5 and p — 3, lasting for about 10^ and 10^ 
periods respectively, while the stickiness time is about 5 x 10^ periods). The 
kind of behavior encountered by these orbits is very sensitive to the choice 
of initial conditions, a fact consistent with the fractal structure of the phase 
space near cantor i. 



34 The Arnold Web of a 4D Map 



Chaotic indicators are widely used in order to visualize the Arnold web of res- 
onances in multidimensional conservative systems. In the following numerical 
examples we consider the 4D symplectic mapping proposed in [20] 



^ _ sin (xj + Vj] 

^"^^ ^ (cos {xj + yj) + cos{zj + tj) + 
T/j+i = Hj + Xj {mod2Ti) 

sin {zj + tj) 

Zj+i - Zj - ^(cos {xj + yj) + cos{zj + tj)+Ay ^ ^ 

tj+i = tj + Zj {mod27T) 

in order to study the behavior of APLE at the chaotic border of a single 
resonance domain. In this mapping, {xj,Zj) are action variables and {yj,tj) 
are their conjugate angles. When e = we have constant values of the actions 

while the values x, z yield also the frequencies, 
i.e., the per step changes of the angles. The lines rriix + m2Z, with mi,m2 
integer, yield the Arnold web of resonances in the action plane. The chaotic 
motions in this mapping along the borders of resonances, when e 7^ 0, were 
studied in detail in [39|20] . Of particular interest are the chaotic motions when 
e is smaller than a threshold value e < eo marking the onset of validity of the 
Nekhoroshev regime [38] (see p9.j). In that case the chaotic motions can be 
identified to the phenomenon of 'Arnold diffusion' [lO] (see [H] for a review). 

Numerical plots of the Arnold web on the action space can be obtained by 
plotting on a color or gray scale the value of a chaotic indicator as a function 
of the initial conditions of the orbits on the action space. Plots of this type, 
for different multidimensional Hamiltonian systems or mappings, were given, 
using different indicators |42f43f30f20f 18|I44] . Here we use the APLE index 
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Fig. 10. (a) A plot of the geography of the resonances in the action space of 
the 4D map ([39]) when e = 0.05, by means of the index Pmax calculated over 
T = 10^ iterations per orbit. The initial conditions are in a 500 uniform grid 
and the initial deviation vector for each orbit is ^ = (■^'O'O'"^)- Each point 
is colored according to the value of Pmax in the range 0.6 < Pmax < 1-4. (b) A 
zoom of (a) around the resonance = 2, produced from a set of 500 x 500 ini- 

tial conditions and the same vector ^ as in (a). The color scale for Pmax is in the 
range 1.5 < Pmax ^ 2. The black points correspond to three initial conditions of 
orbits, namely {xo,yo, zo,to) = (1.67209489949749,0, 0.809297781072027,0) (point 
(1)), (1.67192854271357,0, 0.809236608877722,0) (point (2)), (1.67203944723618,0, 
0.809286279229481,0) (point (3)), (1.67203675,0, 0.809286625,0) (point (4)), 
(1.6720387,0, 0.80928565,0) (point (5)). 

Pmax, for T = 10^ periods, in order to produce a plot of the Arnold web for 
the mapping (139|) . e = 0.05. We set the initial values of the angles yo = to = 
and use a 500 x 500 grid of initial conditions on the action plane (xq, Zq). For all 
these orbits the initial deviation vector is {dxo,dyo,dzo,dtQ) = ("^;0,0,;^). 
The results are shown in figure [TOk . The color scale corresponds to values of 
Pmax in the range 0.6 < Pmax < 1-4, while, if the calculated value oipmax for an 
orbit is outside the above interval, the associated color in the plot is replaced 
by that corresponding to the lower limit 0.6 or the upper limit 1.4. 

The web of resonances is clearly visible in the plot of figure [TOk. which is similar 
to plots of the same system produced in [20] using a different indicator, namely 
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the FLI. Now, these authors studied in detail the diffusion in the action space 
of chaotic orbits starting at the chaotic border of a single resonance domain. 
When e is sufficiently small, the system is in the 'Nekhoroshev regime' in which 
the diffusion coefficient is found to be exponentially small in the inverse of the 
perturbation, i.e., D oc exp ( — (co/e)'') for some positive exponent h (estimated 
as 6 = 0.28 by the same authors). Since the diffusing orbits are in general very 
weakly chaotic, we expect that some of them exhibit the metastable behavior 
associated with a constant rate of production of the Tsallis q-entropy. We 
find this to be the case for some of the orbits explored in [2Dj. If a zoom of 
figure [TUk is made in the neighborhood of a chaotic border of the resonance 
X — 22 = (figure [TOb). the orbits studied in |20j (their figure 5) correspond 
to the initial conditions 

(xo,yo) = (1.67209489949749,0,0.809297781072027,0) (point (1)) or 
(1.67192854271357,0,0.809236608877722,0) (point (2)) and 
(1.67203944723618,0,0.809286279229481,0) (point (3)). The point (3) is near 
a line passing from the center of the chaotic zone of figure [TOb . while the two 
other points are near the edge of the same zone. Figure [TTb then shows the 
time evolution of the APLE, p vs. log^^Q t for these orbits. Clearly, in the case 
of the point (1) in figure [TUb. the resulting orbits shows no visible metastable 
behavior, while such a behavior is clearly exhibited by the two orbits with 
initial conditions on the edge of the chaotic border. In fact, the calculation 
of the Lyapunov exponents A for these orbits [20] shows that the exponent 
of the first orbit is of order A ~ 10~^, while it is two orders of magnitude 
smaller A ~ 10~^ for the orbits near the edge of the chaotic border. We 
conclude that even if a system is in the 'Nekhoroshev regime', not all the 
chaotic orbits manifest the metastable behavior associated with the constant 
rate of production of the q-entropy, but only the orbits with initial conditions 
close to the edge of chaotic borders separating resonance from non-resonance 
domains. 

In order to further substantiate this conclusion, figure [TTb shows the value of 
Pmax for T = 10"^ versus the value of the coordinate xq by a detailed scanning 
along a line passing through the chaotic border of figure [TOb in the direction 
perpendicular to the direction of the resonance under study. We see that as the 
center of the chaotic border is approached the value of Pmax grows abruptly. 
This plot is similar to the plot of Pmax for the crossing of a thin separatrix 
chaotic layer in the 2D standard map (figure [6^), a fact expected since the 
dynamics near the chaotic border of a single resonance domain is qualitatively 
given by a separatrix-like map [28]. At any rate, figure [TTb clearly shows that 
a relatively small value of the APLE (above unity), leading to a long lasting 
metastable behavior of the orbits, can be expected only near the edge of the 
chaotic border of the resonance. In fact, if we take two initial conditions along 
the scanning line of figure [TTb (points (4) and (5) in figure [TOb). the time 
interval of metastable behavior (figure [TTb) is longer for the orbit (4), which is 
closer to the border of the chaotic layer than for the orbit (5) which is closer 
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Fig. 11. (a) The time evolution of the APLE for the three orbits with initial 
conditions as in figure [TOb (thick dotted black curve for the point (1), gray curve 
for the point (2), and thin black for the point (3) point), (b) The value of Pmax 
after T = 10^ iterations as a function of the initial condition xq of a set of orbits 
in a segment crossing perpendicularly the resonance line of figure [TOk . (c) Time 
evolution of p for two orbits in the same line, with initial conditions as in figure 
llOb . The thick dotted line corresponds to the point (4) (orbit near the border of the 
resonance chaotic layer) and the thin continuous line to the point (5) (orbit near 
the center of the resonance chaotic layer). 
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Fig. 12. As in [TOb. but for the average value < p > calculated over 10"^ iterations 
per orbit. 

to the center of the layer. However, both orbits are relatively closer to the 
center than the orbits (1) and (2), and, consequently, they both have shorter 
metastable time intervals than the latter orbits (figure [Tib) . 



Finally, figure [I2] shows the average value < p > over an interval < t < 10^ 
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for all the orbits with initial conditions in the same interval as in figure [TTb . 
These orbits have a value of < p > always above unity, and the mean of all 
values in the interval is < p > = 1.27, yielding an entropic index q = 0.21. 



4 Conclusions 

We studied the connection between the production of Tsallis q-entropy and 
the behavior of the variational equations of motion for weakly chaotic orbits 
in conservative dynamical systems. Our main findings can be summarized as 
follows: 

1) The solutions of the variational equations present long transient time inter- 
vals in which the length of the deviation vector increases almost as a power-law. 
This allows to define an almost constant average power-law exponent (APLE) 
during the whole transient interval. 

2) This 'metastable' behavior can be justified theoretically by showing that 
it is caused by the growth of the deviation vectors inside separatrix-like thin 
chaotic layers, which is of the form ^(t) ^ at + e'^*, with a » 1 and A << 1. 
The latter law appears almost as a power law ^(t) ^ for time intervals up 

to t2c ~ 

3) The average value of the APLE in a thin chaotic layer can be used to de- 
termine an average value of the q-entropic index for which the Tsallis entropy 
exhibits a constant rate of increase. 

4) The APLE can be used as an efficient numerical indicator distinguishing 
regular from weakly chaotic orbits. In that respect, it is equally powerful to 
other established indicators, as the FLI or the MEGNO. The advantage of the 
APLE is that it gives also the average value of the q-entropic index within a 
weak chaotic layer. 

5) Numerical implementations of the APLE are given in low-dimensional sym- 
plectic mappings. The APLE is calculated in a thin chaotic layer and in the 
stickiness region of an island of stability in the 2D standard map. We then use 
it in order to visualize the Arnold web of resonances in a 4D map, and calcu- 
late the q-entropic index at the chaotic border of a single resonance domain 
of the same map. 
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